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We present the results of a combined experimental and theoretical investigation of the 
surface-tension-driven coalescence of flexible structures. Specifically, we consider the dy- 
namics of the rise of a wetting liquid between flexible sheets that are clamped at their 
upper ends. As the elasticity of the sheets is progressively increased, we observe a system- 
atic deviation from the classical diffusive- like behaviour: the time to reach equilibrium 
increases dramatically and the departure from classical rise occurs sooner, trends that we 
elucidate via scaling analyses. Three distinct temporal regimes are identified and subse- 
quently explored by developing a theoretical model based on lubrication theory and the 
linear theory of plates. The resulting free-boundary problem is solved numerically and 
good agreement is obtained with experiments. 
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1. Introduction 

Surface tension plays a role in an enormous variety of industrial problems and natural 
phenomena: the string of dew drops on a spider web, the spontaneous imbibition of water 
into soil, the stability of water-walking arthropods, the wet lithography of micro-electro- 
mechanical systems, the rise of sap in trees, and the wet adhesion of insects on plants, to 
name just a few. It is clear that the elasticity of the substrate is relevant to understand 
some of these phenomena, including the clumping of the flexible bristles of tarsi of insects 
when tarsal oil is released (Eisner & Aneshansley 2000 ) or the adhesive failure (stiction) of 
micro-cantilevers bent by capillary forces at the evaporating menisci ( Mastrangelo fc Hsu| 



1993). Although static problems involving the interplay between capillarity, elasticity. 



and possibly other forces are well studied, there are very few studies of the dynamics of 
elastocapillarity, which is the subject of this paper. 

Early studies of elastocapillarity dealt with the deformation of a solid in contact with a 
liquid drop. From the viewpoint of mechanics, when a liquid drop adheres to a substrate, 
there must be a balance of the forces acting at the contact line. The familiar Young's 



equation can be interpreted as a balance of the tangential forces. Some authors (Lester 
1961 Fortes|p!984 ) suggested that the normal component is balanced by a force exerted 



by the solid, which deforms upon contact. This deformation is generally negligible on a 



rigid substrate, but can lead to the formation of a ridge (Andrade et al. 1979 Pericet- 



Camara et al. 2008) or wrinkles (Huang et al. 2007) on a soft substrate. Motivated by 



the recent interest in soft materials in engineering, nano- fabrication, and biomechanics, 
several authors have addressed the question of how slender objects deform under capillary 



forces dBico et a/.||2004| |Kim fc Mahadevan] [20061 |Py et a/.|[2QQ7l [Pokroy et a/.||2009D , 
mostly considering equilibrium configurations. 
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Nevertheless, some elastocapillary effects are essentially dynamic (e.g. sap flow, spore 
release), and it is of fundamental and practical importance to estimate the time to reach 
equilibrium in such systems, as well as to understand their transient dynamics. Recent 
examples of such elastocapillary dynamics include studies of capillary-driven flow in de- 
formable porous media ( Siddique et a/._2QQ9) and in fle xible channels where gravitational 
effects could be neglected ( van Honschoten et a/.|2QQ7 Aristoff et a/.|2QlQ ). Here, we con- 
sider a model system to study elastocapillary dynamics with gravity, namely capillary 
rise between flexible sheets that are clamped at their upper ends. In doing so, we extend 
the work of Bico et al. ( 2004 ) and Kim & Mahadevan ( 2006 ) by considering the dynamics 
of elastocapillary rise. The experimental set-up and results are described in ^ We begin 
with statics, identify three different regimes, and then address the transient behaviour. 
We show how the evolution of the position of the meniscus departs from that of classical 
capillary rise (i.e. rise between rigid boundaries) and exhibits an unusual behaviour: the 
meniscus slows down, reaches a plateau, then reaccelerates before finally stopping. In ^ 
we formulate a free-boundary problem to model the dynamics, and in Q we describe 
the solution technique and present the numerical solutions. The transient dynamics is 
analysed in ^in light of experimental and numerical results. 



2. Experiments 

2.1. Experimental setup and dimensionless numbers 

In this section we describe our experimental setup to investigate the dynamics of capillary 
rise between flexible sheets (Fig. [T]). Two vertical glass sheets (length ^, thickness h = 
160 /im, width w = 5 mm and bending stiffness per unit width B = (2.2 ± 0.2) x 10~^ 
N-m) are initially parallel. The sheets are clamped at their upper ends, with a separation 
distance 2ho, and free at their lower ends. At time t = 0, the lower ends of the sheets are 
brought into contact with a bath of silicon oil (viscosities and densities /i = 0.096 Pa-s , 
p = 963 kg-m-3 (VlOO), fi = 0.048 Pa-s, p = 960 kg-m-^ (V50), p = 0.0093 Pa-s, p = 
930 kg-m~^ (VIO), and surface tension 7 = 0.021 N-m~^) that perfectly wets the glass 
(i.e. the equilibrium contact angle Oe = 0). The deflection of the sheets and the position 
of the meniscus as it rises between the sheets are recorded from the side with a digital 
camera and tracked using customised image- analysis software in MATLAB®. 

We take the length of the sheet i and half the initial gap ho as, respectively, the 
characteristic lengths for the meniscus position Zm{t) and the deflection of the sheets 
from the vertical h{z^t). We can then estimate the capillary energy Ec ~ w£^, the 
gravitational energy Eg ~ pgwh^f' and the bending energy E^ ^ Bwh^/i^ for our 
system, where g is the gravitational acceleration and here we use ~ to indicate order-of- 
magnitude scaling estimates for which prefactors are typically order one. Balancing Ec 
with Eg^ Ec with Ee^ and Eg with Ec yields three characteristic lengths: 

4. = ^, 4c=f^)''',and ieg=(^Y\ (2.1) 

The capillary-gravity length leg is the classical equilibrium height for the capillary rise 
between two rigid plates. The elastocapillary length l^c corresponds to the minimal length 
above which capillary forces can bring the sheets together and has been identified previ- 



ously in studies of the failure of micro- mechanical structures upon drying ( Mastrangelo 
fc Hsu|1993||de Boer fc Michalske|19"99l ) or hair clumping ^Bico et a/.|2004| ). The scale I eg 
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Figure 1. Sketch of the experimental setup. The top hatched region corresponds to the clamped 
edge of the glass sheets. The shape of the meniscus is indicative of a typical experimental 
observation. We track the position of the apex of the meniscus. 



is an elastogravity length. Our system is thus described by two dimensionless numbers, 

I 



S=-^and^=f-^) , (2.2) 



^cg \-^. 

where B is the Bond number and is the elastocapillary number. Note that we may also 
define an elastogravity number Q = {£/£eg)^ = Bco . By changing the spacing 0.18 mm ^ 
^ 0.39 mm, we varied the elastocapillary length 13.4 mm ^ l^c ^ 20 mm. The 
elastocapillary number S' is varied within a wide range (0.18 ^ ^ 192) by changing 
the length of the sheets 10 mm ^ ^ ^ 50 mm. The value of the Bond number B is also 
affected by the changes in ^, but remains of order one (1 ^ S ^ 7) for all experiments. 
We note that for S ^ 1, the meniscus will fill the entire gap between the sheets (i.e. 

2.2. Equilibrium configuration 

When the free ends of the sheets are brought into contact with the bath, the liquid spon- 
taneously rises up, and the induced pressure distribution leads to an inward deflection 
of the sheets until the meniscus stops at a finite height. This equilibrium state depends 
on the combined effects of surface tension, gravity, and elasticity. Depending on the two 
dimensionless parameters (^, S), three different configurations are observed, as depicted 
in Fig. |2j For (o < 10, the sheets are slightly deflected while the lower ends remain open 
(Regime I). For > 100, the sheets deflect and coalesce over a distance comparable to 
their length (Regime III). Nevertheless, the sheets do not make dry contact; a thin liq- 
uid film remains between them. The statics of these two regimes have been investigated 



previously (Bico et al] 20041 iKim & Mahadevan 2006). In addition, we have identified 



an intermediate regime (Regime II) wherein the sheets touch at the tip, but there is a 
nonzero angle between them. 

The equilibrium shape of the sheets is found by solving numerically the nonlinear free- 



boundary problem developed by Kim & Mahadevan (2006), with a different boundary 
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Figure 2. Equilibrium configurations. Regime diagram of elastocapillary rise in the parameter 
space (cf, B) obtained by studying the static problem experimentally (symbols) and numerically 
(shaded regions). Three regimes, labeled I (o), II (o) and III (□) are identified and illustrated 
with experimental photographs and numerically deduced profiles for (I, •) S' — 7.2, B — 1.7, 
(II, ♦) ^ = 12, B^ 1.9 and (III, ■) ^ = 192, B = 3.9. In the region labeled I/II, we predict 
numerically the coexistence of two regimes (I and II), in agreement with our experimental 
observations. 



condition at the lower ends of the sheets necessary to find an equiHbrium shape for 
intermediate S' ^ which corresponds to Regime II (see ^and Table [T]). We compiled the 
permissible equilibrium configuration (s) for 120,000 different pairings of S' and S, and 
present the results in Fig. [2] (shaded regions). Good agreement with our experimental 
observations is observed. We note that for 10 < < 30, solutions for regimes I and II 
coexist. 



Kim & Mahadevan ( 2006 ) performed a scaling analysis to determine the onset of large 



deflection of the sheets, which can be interpreted as the upper boundary of Regime I. To 
do so they balanced the torque (per unit width) exerted on the sheets, Bho/i'^^ with the 
capillary-gravity torque, pg£lg{£ — £cg/'^)' They considered the limit i ^ icg^ (i-e. B 1), 
wherein the capillary-gravity torque becomes pgi^gi^ and thus they obtained a boundary 
that is given hy S' ^ B. Our numerical results suggest that this scaling is indeed valid for 
large Bond numbers (not shown). However, our experiments and numerically predicted 
equilibrium shapes show a different behaviour for B < 10, which can be understood 
by considering the limit i/icg 1 (i-e. S ^ 1), wherein the capillary-gravity torque 
becomes pgi^. Hence, we find a boundary that is given by ~ in agreement with 
our numerical solutions shown in Fig. |2] 

2.3. Dynamics 

We now focus on dynamics, in particular the time to reach these equilibrium config- 
urations. We measured the time t* taken by the meniscus to reach 99% of its final 
height. Results obtained for ho = 0.18 mm and three different viscosities are presented 
in Fig. [3]^ a), where we varied the length of the sheet i and hence the elastocapillary 
number (o. For short sheets, is roughly constant and on the order of a minute. The 
sheets deflect slightly and the final state is close to the one obtained for rigid boundaries 
(Regime I). For larger sheets, the deflection becomes important (Regimes II and HI) and 
t* increases rapidly with increasing i. 

To estimate the time to reach equilibrium, we assume that the speed of the meniscus 
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Figure 3. Time to reach equilibrium, (a) Time t* in seconds (s) required to reach 99% of the 
equiUbrium height versus the length of the sheet i in millimeters (mm). Data points correspond 
to experiments with a constant initial gap /lo = 0.18 mm and with silicone oil of viscosity VIGO 
(♦), V50 (*) and VIO (+). (b) Time t* scaled by the characteristic time to reach equilibrium in 
the classical rigid case, Tr, versus B^<§ . The dotted line corresponds to t* = 3.6Tr and the solid 
line to t*/rr = ^ ^ i.e. t* = r^e- 



u is given by that for a Poiseuille flow u = /S.ph^/ {^tfit). Provided the deflection of 
the sheets is small, the pressure difference across the meniscus may be approximated 
by Ap = and the characteristic length is the capillary-gravity length leg- The 

time to reach equilibrium is thus given by a balance be tween capillary, viscosity, and 



1997), which is independent of the 



gravitational forces, leading to = ^3^^^^^ 2 (|Quere 
length of the sheets I. 

When the deflection of the sheets is appreciable, the two dominant effects are viscosity 
and elasticity. The pressure difference is set by the typical deflection of the sheets h^^ 
i.e. Bho/i^ = Ap. Therefore, we can deflne a characteristic visco-elastic time r^e = f = 

= S^^fr^, which rapidly increases with increasing £. 

The time to reach equilibrium, rescaled as t*/rrj is presented in Fig.pfb) versus B^(o. 
The data follow two different behaviours. For B'^S' < 10, we recover the classical result for 
capillary rise between rigid sheets: t* ^ indicating that elastic effects are negligible. 
Conversely, for B^(o ^ 10 the time to reach equilibrium scales with the visco-elastic time, 
i.e. t*/r^ = S^^f , indicating that elastic effects are dominant. 

Next, we examine the evolution of the shape of the sheets and the position of the 
meniscus for three experiments representative of Regimes I, II and III (Fig. [4|. The 
dynamics of imbibition between flexible sheets is compared to the classical case of rise 
between rigid sheets (black lines in Fig. [4|, wherein the meniscus f ollows a diffusive-like 



behaviour at early times, i.e. cx (Bell & Cameron 1906), and is then slowed 



by gravity until eventually reaching its equilibrium height leg- Iii our experiments (grey 
points in Fig.[4|, the meniscus follows this classical behaviour at early times, as the sheets 
do not appreciably deflect. At later times, the sheets begin to bend and the meniscus 
position deviates from that of classical capillary rise. For small (o (Fig. [4]^ a)), the sheets 
are slightly deflected. The meniscus slows down and stops at a higher position than 
the classical height. For intermediate (o (Fig. |4jb)), the ends of the sheets touch during 
imbibition. The meniscus slows down, reaches a plateau, then surprisingly, accelerates 
until flnally stopping. For large (o (Fig. [4]^c)), the sheets touch soon after imbibition 
begins, when the meniscus is quite close to the bath and gravity has not yet contributed 
to its deceleration. At intermediate times, the position of the meniscus dips below that 
of classical rise and the meniscus slows down considerably. After roughly one minute 
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Figure 5. Rescaled position of the meniscus (zrn/^)(2<f )~^/^ as a function of time 
T — t/Tve — t Bhl/{3iJ.i^). The sohd Une corresponds to a diffusive-Uke behaviour Zm oc \fWf . 
Inset: evolution of the cross-over time T' at which the meniscus position deviates from the clas- 
sical prediction as a function of for experiments (o) and numerics (+). Numerical solutions 
are described in ^ The solid line represents the predicted scaling T' oc . 



has elapsed, the sheets begin to coalesce and the meniscus accelerates until eventually 
reaching its equilibrium height. 

To obtain a better understanding of the dynamics of elastocapillary rise, we compare 
the motion of the meniscus for several experiments as the elastocapillary number S is 
progressively increased. Since the meniscus position at early times is given by oc t^l^ ^ 
or equivalently, in dimensionless terms, Zm = ^m/^ oc ^/^(t/r^e)^^^, we use the char- 
acteristic length and time scales M'^^'^ and r^g, respectively, to rescale the experimental 
data. The evolution of the dimensionless position of the meniscus (z^/-^)(2^f )~^/^ with 
time T = t/Tye is presented in Fig.js] The family of curves obtained for various (o values 
collapse onto the same line at early times, corresponding to a diffusive- like behaviour of 
the advancing meniscus, i.e. Zm oc \/ ST . Here, we have used ex to distinguish order-of- 
magnitude results from results indicating proportionality. We have also verified that the 
fluid viscosity enters the dynamics only through r^g, as shown by the collapse, in Fig.jS) 
of the curves for three different viscosities (VIO, V50, and VlOO) at = 11.4. 

The position of the meniscus deviates from T^/^ at a time T' . Increasing S decreases 
this cross-over time T', as the sheets bend more easily. Neglecting the deceleration due to 
gravity, we expect this deviation to occur when the deflection of the sheets due to capillary 
forces becomes important, that is for Zm ^ ^ec-, oi" in dimensionless terms, Zm ^ S~^l^ . 
Before reaching this value, the meniscus evolves accordingly to T cx Z^j S . Hence, we 
expect a cross-over time of T' oc . The dependence of T' on S is given in the inset 

of Fig. [5] The predicted scaling is indeed recovered. For times greater than the cross-over 
time, the speed of the meniscus slows down dramatically. As the elastocapillary number 
S is increased, the motion of the meniscus becomes increasingly complex. In particular. 
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we observe a plateau region followed by an acceleration of the meniscus. In order to 
understand this behaviour quantitatively, we now formulate a theoretical model for the 
dynamics of elastocapillary rise. 



3. Theoretical formulation 

In this section we formulate a free-boundary problem that describes the evolution 
of the shape of the sheets and the rise of liquid between them. We first consider the 
deflection of the initially vertical sheets. Owing to the symmetric conflguration, we need 
only consider one of the sheets. Provided that the initial gap, 2/io (set by the clamping 
distance), is much less than the length of the sheet, ^, the slope of the sheet, h^/l^ is 
small. We assume that the sheet is sufficiently thin, 6 <C <C ^, so that we may use the 
linear theory of beams to describe the deffection, h{z^t)^ from the vertical. Hence, 

Bhzzzz + rnhtt = p{z, t), (3.1) 

where subscripts denote partial derivatives, B = Eb^ /12{1 — f'^) is the bending stiffness 
per unit width, E the Young's modulus, v the Poisson ratio, m the mass per unit length 
per unit width, and p(z, t) the force per unit area (pressure) acting on the sheet. The 



second term in (3.1 ) arises from the inertia of the beam, and may be neglected if the time 



scale of the flow, r^g? is much longer than the reaction time of the beam to an applied 
load, r = \J obtained by balancing the characteristic bending and inertial stresses. 
We use this quasi-static description of the sheet wherein r^g ^ Hence, in the following, 

Bh,,,,=p{z,t). (3.2) 

We note that the weight of the sheet gives rise to a stretching strain of order mg£/Eb. 
Comparing this to a typical bending strain Bho/i'^ leads to a maximum length ig ^ 

1 /3 

(^Eb'^ho/mg) , below which stretching effects can be ignored. In our experimental study, 
we have £ <^ is — cm. 

We now consider the incompressible, viscous flow of a perfectly wetting liquid {6e = 0) 
that enters the region between the sheets owing to a reduction in pressure at the meniscus 
relative to the bath. Our results are extendable to partially wetting liquids (0 < < 
7r/2). In what follows, we take the limit ho <C relevant to our experiments, so that 
edge effects may be safely neglected. Hence, the pressure drop is given by where 

= cos 6d/h{zm,t) is the curvature at the meniscus z = Zm{t) and 6d is the apparent 
dynamic contact angle. Dynamical features of the moving contact line are treated by 
using a semi-empirical relationship for the case of total wetting, in which the dynamic 
contact angle, Od^ is deflned as 



as flrst measured by Hoffman (1975) and later established by Tanner ( |1979 ). This rela- 



tionship accounts for viscous dissipation at the contact line and leads to a reduced driving 
pressure during imbibition ( Guyon et a/.||2QQT ). The parameter F = 6\n{xM /xm) arises 



from evaluating the viscous dissipation near the contact line, where it is necessary to 
introduce the upper xm and lower Xm cut-off lengths to prevent the divergence of the 
integral measuring the total viscous dissipation. No exact solution for the microscopic 
description is currently available, and thus we treat F as an adjustable parameter. In 
order to determine its value for our experimental study, we performed a series of control 
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experiments where we examined the imbibition of sihcone oil between rigid glass bound- 
aries, both in planar and axisymmetric configurations, and with and without gravity. 
The value of T which best describes the results of these control experiments is F = 60, 
and so we take it to have this value in our numerical simulations of elastocapillary rise. 
This value of F is comparable with those found in the literature. Finally, we note that 
our numerical results are only slightly dependent on the value of F; simulations with 
50 < F < 70 give similar results. 

A one-dimensional approximation to the viscously dominated fiow is given by Darcy's 
law, 

dp 

where u is the cross-sectionally averaged liquid velocity, p the local pressure and z the 
axial coordinate. The permeability, k = /i^/3, where 2h is the total distance between 
the sheets, is based on the lubrication approximation ( Guyon et a/. [2001, ) which requires 
\hz\ <C 1, and is consistent with our description of the sheets. 



Differentiating (3.2) once in z and combining the result with (3.4) gives an expression 
for the liquid velocity, 

u{z, t) = -j^ {Bh^zzzz + pg) , (3.5) 

which may be evaluated at z = Zjn(t) to find an expression for the speed of the meniscus. 
Next, combining (3.5) with the statement of mass conservation, ht + {hu)z = 0, yields a 



nonlinear evolution equation for the defiection of the sheets in the liquid-filled region: 

+ 3M.)- (3.6) 



ht = -- {Bhh 
3/i 



zzzzzz 



Let h{z^t) denote the defiection of the sheets in the liquid- free region. Using (3.2) with 
p{z^ t) = 0, we have 

Bh,,,, = 0, (3.7) 



which may be integrated analytically. Equation (3.6) is similar to the equation studied 



by Hosoi & Mahadevan (2004), who considered the dynamics of an elastic sheet that was 



clamped at one end and lubricated from below by a single fluid. Our study is distinguished 
by the influence of surface tension and the presence of two fluids whose boundary advances 



according to (3.5) and depends on the local shape of the sheet. 



The evolution equation and boundary conditions may be rendered dimensionless using 
the characteristic length and time scales presented in ^2.1 Letting H = h/ho^ Z 
Zm = Zm/^j and T = t/rye^ transforms (|3.5[), evaluated at the meniscus, to 



z/e, 



dZm 



zzzzz 



SB) 



Z — Zrr 



(3.; 



and (3.6) to 



Ht = H^Hzzzzzz + SH^HzHzzzzz + ScoBH^Hz- (3-9) 
Two initial conditions and six boundary conditions are necessary to solve the free- 



boundary problem specifled by (3.8) and (3.9). The flrst two boundary conditions are 
Hzz{0,T) = 0, and Hzzzz{^-,T) = 0. These correspond, respectively, to zero moment 
and zero pressure (relative to the atmosphere) at the lower end of the sheets. If the 
free ends of the sheets remain open (Regime I), we have the condition i^zzz(0,T) = 0, 
which corresponds to zero force at the lower end of the sheets. If the sheets are in 
contact (Regime II), they exert an equal but opposite force on each other, and the 



10 



C. Duprat, J. M. Aristoff and H. A. Stone 



Table 1. Boundary conditions for (3.9) in each of the three temporal regimes: Regime I (open 



ends), Regime II (ends in contact), and Regime III (ends coalescing). 



1 Regime I 


Regime II | 


Regime III | 




Hzz{0,T) = 


-- 1 




Hzzzz{0, T) 


= 1 


\Hzzz{0,T) = Q\ 1 1 




HiO,T) = 0\ 


if(Z^Ze,T) = 1 


1 1 1 Hz{Z ^Zc,T)=0 1 


1 Hzzzz{Zm,T) = -£" cosed/ H{Z^) \ 


\H{Zm,T) = l- 


^{HzZz{Zm,T)){Zm — 


l)' + |i/zz(^m,T)(Z^-l)2| 


1 Hz{Zm,T) = - 


-^{HzZz{Zm,T)){Zm - 


-lf + Hzz{Zm,T){Zm-l) 1 





force- free condition i^zzz(0,T) = must be replaced by i^(0,T) = 0. If the sheets 
are in contact over a finite length < Z ^ ZciT)^ where Zc = ^c/^ and Zc is the 
highest point at which the sheets contact, as in Regime III, we have an additional un- 
known, Zc, and an additional condition Hz{Zc^T) = 0, that enforces a tangentially 
smooth contact between the sheets. At the meniscus, the jump in pressure given by 



(3.2) is equal to the Laplace pressure 7/^:, which can be written in dimensionless form as 
HzzzziZm^T) - HzzzziZm.T) = COS ed/H{Zm). 

The remaining two boundary conditions for H{Z,T) are obtained by expressing the 
matching conditions at the meniscus between the liquid-filled and the liquid-free regions 
in terms of the boundary conditions at the top end of the sheet (fixed and clamped, 
i.e. H{1,T) = 1 and Hz{l^T) = 0). Combining the conditions of equal deflection 
^(Z^,T) =i7(Z^,T), equal slope ^z(^m,T) =i7z(Z^,T), equal moment fe(Z^,T) 
Hzz{Zm-,T) and equal shear force (since the liquid is perfectly wetting) Hzzz{Zm^ T) = 
Hzzz{Zm,T), yields 

H{Zm,T) = 1 - ^{Hzzz{Zm:T)){Zm - if + ^fc(^m,T)(Z^ - 1)^ (3.10) 

and 



HziZm.T) = --{HzZz{Zm,T)){Zm - if ^ H z z{Zm,T){Zm - 1). (3.11) 



In doing so, the solution of (3.9) for the deflection of the sheets in the liquid- filled region 
automatically satisfies the boundary conditions at Z = 1, and prescribes the deflection of 
the sheets in the liquid-free region Zm < Z ^ 1. The boundary conditions are summarized 
in Table [l] Finally, we take for the initial conditions H{Z^O) = I and Z^(0) = 0.001. 
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4. Solution of the free-boundary problem 

The solution to the free-boundary problem defined in ^is divided into two stages. For 
early times (T ^ 1) the defiection of the sheets is small, time enters only through the 
boundary conditions at Z^(T), and the free-boundary problem admits the approximate 
solution 

"^'■^^ ^ ' - iS; + (f - f + 1) (t - § - s) 

for the defiection of the sheets in the liquid-filled region, and the liquid advances as though 
it were between rigid sheets. Hence, with Zm = ZmB = Zm/^cg and T = T(^B^ = t/r^, 
we find 

dZm _ COS Od _ ^ ,^2) 
df ~ Zrr,, ' 



in accordance with classical capillary rise. When cos^^i = 1, (4.2) can be integrated to 
give an implicit expression for Z^(T): 

T = -Z^-ln(l-Z^). (4.3) 



At early times, (|4^ reduces to Z^(T) = V2T, or equivalently, Z^(T) = VWT for 
Zm <^B~^. When cosOd 7^ 1, the meniscus is slowed down at early time and reaches the 
classical behaviour at late times. 

For late times, when the defiection of the sheets is significant, we solve numerically 



the free-boundary problem using a method similar to the one described in |Aristoff et al 



(2010). We transform the time-dependent domain ^ Z ^ Zm{T) to the fixed do- 



main ^ 5* ^ 1 via the substitution S = Z/Zm{T)^ and rewrite the time and spatial 



derivatives in (3.8) and (3.9) accordingly. Then, the resulting equations and boundary 
conditions are discretized using an implicit finite-difference scheme and solved numeri- 
cally in MATLAB®. We use for initial conditions the profile predicted by (4.1) and the 



meniscus speed given by (4.2). 



The evolution of the position of the meniscus obtained from our numerical simulation 
is shown in Fig. [6j Comparison with experiments (Fig. [5| shows that our model cap- 
tures the essential features of the experiments, the collapse onto \^WT at early times, 
and the observed plateau and subsequent acceleration at late times and for large (o. 
At early times, the position of meniscus follows the expected diffusive-like behaviour 
Zm = aV2lT. The pre-factor a = 0.65 < 1 arises solely from the deceleration induced 



by the dynamic contact angle (3.3) and is identical in the experiments and in the numer- 
ical solutions. Direct comparison between the experimental and numerical results (inset 
in Fig. [6| shows good quantitative agreement between the experimental data and the 
corresponding numerical solutions. The discrepancy at late times may result from the 
difficulty in accurately measuring the bending stiffness B or the presence, in Regimes II 
and III, of a thin liquid layer between the sheets, which is neglected in the numerics and 
effectively increases the gap in the experiments. We also verify the scaling law for the 
cross-over time T' (inset in Fig. |5| with a numerical pre-factor of 0.25 corresponding to 
a cross-over when the meniscus reaches Zm — 0.5 iec Foi" > 35, the equilibrium shape 
is not achieved in the simulations, owing to the meniscus speed reaching nearly zero in 
the plateau region. 
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Figure 6. Rescaled position of the meniscus (zrn/^)(2^f ) ^ as a function of time t/rv 

for numerics. The sohd hne corresponds to — a\j2<§T with a pre- factor of a 
the dynamic contact angle effect given by (3.3). Inset: direct comparison between experimental 
data (grey points) and numerical solutions (black lines). 



0.65 due to 



5. Discussion 

The interpretation of our results may be made by considering the geometry of the 
channel. As the sheets are deflected, the gap evolves from being uniform to having a 
diverging geometry when the sheets touch. The asymptotic behaviours of the position 
of the meniscus are shown in Fig. [7^ (experiment) and Fig. [7]3 (numerics) for matching 
values of S and B. We first consider the experimentally observed dynamics (Fig.[7]^a)). At 
early times, the meniscus obeys the diffusive law oc T^l^ expected for a uniform gap 
in the absence of gravity. When T = T' ^ the deflection of the sheets becomes important 
and the gap resembles a wedge. The meniscus follows the power law expected in a wedge 
geometry Zm cx T^/^ ( Warren|2QQ4 ). We then observe a plateau of duration Twait during 



which the meniscus advances slowly. The channel formed by the two sheets adopts a 
complex geometry, and the data roughly follows a power law oc T^/^^. This behaviour 
is consistent with imbibition in a diverging channel having shape h{z) = h p + 0z^, where 



the time-dependence of the meniscus position is given by Zm oc 7^1/(2^+1) ( Reyssat et al 
|2QQ8| ). Our results are close to n = 6, although any value of n > 1, hence an exponent of 
l/(2n + 1) < 1/3, is permissible. We note, however, that gravity will slow down the rise 
of liquid and so the analogy is qualitative at best. 

The asymptotic behaviours found experimentally are recovered in the numerical results 
(Fig. [rj^b)), which allow us to accurately investigate the changes in the shape of the 
sheets during imbibition. In particular, we identify the time at which the sheets touch 
(T = 1.18- 10~^) as the time at which the meniscus reaches a plateau, before accelerating 
owing to the new geometry of the channel. We obtain good agreement between the 
experimental data and the corresponding numerical solution. Nevertheless, the duration 



Dynamics of elastocapillary rise 



13 




10 



10"' 10' 10° 



10' 




10" 

Figure 7. Position of the meniscus (zm/'^)(2^f )~"^/^ versus time for = 11.4, 23 = 2 for (a) 
experiments and (b) numerics. The sohd hues correspond to Zm = 0.65^2^ T^/^ When T = T\ 
the meniscus follows the power law expected in a wedge geometry Zm = 0.28a/2^ T^^^ (dashed 
lines). At later times, the channel formed by the two sheets adopts a complex geometry and the 
meniscus roughly follows a power law Zm oc T^^^^ (dash-dotted lines) with a pre- factor of 0.12 
(experiments) and 0.11 (numerics). Insets show (a) processed experimental photographs and (b) 
computed profiles. The corresponding times are denoted by the symbols (•) on the grey curves. 



of the plateau Twait is overestimated in the numerical version of the imbibition, due to 
the sensitivity of the numerical scheme as the speed of the meniscus approaches zero. 

In closing, we have provided a framework to study the dynamics of capillary-driven 
flow between flexible boundaries, with particular attention given to a model system where 
capillarity, elasticity, and gravity compete. We have identified the relevant length and 
time scales, and shown that the time to reach equilibrium sharply increases with increas- 
ing elasticity of the boundaries. Since so many natural phenomena involve elastocapillary 
effects, it is reasonable to think that the dynamical features we describe here will be use- 
ful for rationalizing the time-dependent behaviour of a variety of biological and industrial 
systems. 
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